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ABSTRACT 

We complete our study of pulsars' non-uniform surface temperature and of its effects 
on their soft X-ray thermal emission. Our previous work had shown that, due to the effect 
of gravitational lensing, dipolar fields cannot reproduce the strong pulsations observed 
in the four nearby pulsars for which surface thermal radiation has been detected: PSR 
0833-45 (Vela), PSR 0656+14, PSR 0630+178 (Geminga), and PSR 1055-52. Assuming 
a standard neutron star mass of 1.4 M & , we show here that the inclusion of a quadrupolar 
component, if it is suitably oriented, is sufficient to increase substantially the pulsed 
fraction, Pf, up to, or above, the observed values if the stellar radius is 13 km or even 
10 km. For models with a radius of 7 km the maximum pulsed fraction obtainable, 
with (isotropic) blackbody emission, is of the order of 15% for orthogonal rotators (Vela, 
Geminga and PSR 1055-52) and only 5% for an inclined rotator as PSR 0656+14. Given 
the observed values, this may indicate that the neutron stars in Geminga and PSR 
0656+14 have radii significantly larger than 7 km and, given that very specific quadrupole 
components are required to increase Pf, even radii of the order of 10 km may be unlikely 
in all four cases. However, effects not included in our study may possibly seriously 
invalidate this temptative conclusion. 

We confirm our previous finding that the pulsed fraction always increases with photon 
energy, below about 1 keV, when blackbody emission is used and we show that it is due to 
the hardenning of the blackbody spectrum with increasing temperature. The observed 
decrease of pulsed fraction may thus suggest that the emitted spectrum softens with 
increasing temperature and that this observed effect must be of atmospheric origin. 

Finally, we apply our model to reassess the magnetic field effect on the outer boundary 
condition used in neutron star cooling models and show that, in contradistinction to 
several previous claims, it is very small and most probably results in a slight reduction 
of the heat flow through the envelope. 

Submitted to The Astrophysical Journal. 

Subject headings: stars: neutron pulsars: individual: Vela, 0656+14, Geminga, 

1055-52 X-rays: stars dense matter stars: magnetic fields 
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1. INTRODUCTION 

In a previous paper (Page 1995a; hereafter Paper 
I) we had presented a simple but realistic model of 
temperature distribution at the surface of a magne- 
tized neutron star. This model was used to study 
the effects of temperature inhomogeneity on the neu- 
tron star thermal emission and compare them with 
the observed properties of the four neutron stars for 
which this emission has been detected by the ROSAT 
satelite (Ogelman 1995a; Page 1995b): PSR 0833-45 
(Vela; Ogelman, Finley, & Zimmermann 1993), PSR 
0656+14 (Finley, Ogelman, & Kiziloglu 1992), PSR 
0630+178 (Geminga; Halpern & Holt 1992) and PSR 
1055-52 (Ogelman & Finley 1993). As a first step 
we did not take into account magnetic effects in the 
atmosphere which are also very large. We also re- 
stricted ourselves to dipolar magnetic fields and the 
three main results we obtained were: 

— Magnetic effects on heat transport in the neutron 
star envelope do induce very large temperature dif- 
ferences at the surface but when gravitational lens- 
ing is taken into account and only dipolar fields are 
considered the predicted pulsed fractions are smaller 
than those observed in the soft band where the surface 
thermal emission is seen. In particular, in the case of 
1.4 Mq neutron stars with small radii ~ 7.0 — 9.0 km, 
the pulsed fractions Pfs are below 1%. With very 
small radii < 6.5 km, and the same mass, the gravi- 
tational beaming can however increase Pf up to 6 - 
7%. Observed values of Pf range between 10 - 30%. 

— With dipolar fields the observable light curves 
are very symmetrical and their shapes do not corre- 
spond to the observations. Together with the previ- 
ous result, this lead us to conclude that the inclusion 
of dipolar fields only is not sufficiently adequate to 
model the surface magnetic field of the four observed 
neutron stars. 

— The amplitude of the pulse profile increases with 
increasing photon energy. This result does not corre- 
spond to what is observed in the soft X-ray band in 
the cases of PSR 0656+14, Geminga and PSR 1055- 
52. 

We complete here this study and consider also the 
effects of these surface temperature distributions on 
the cooling of neutron stars. Since the completion of 
the work presented in Paper I, complementary results 
considering the magnetic effects in the atmosphere, 
where the emerging spectrum is generated, have been 
presented (Pavlov et al. 1994; Zavlin et al. 1995) 



and clearly showed that they can be as important as 
the magnetic effects in the envelope that we consider 
here and in Paper I. A complete study must obviously 
include both aspects of the problem, envelope and at- 
mosphere. However, as long as the exact chemical 
composition of a pulsar's surface is not known, anal- 
ysis with blackbody (BB) emission will still remain a 
mandatory first step and a reference to which atmo- 
sphere models will be compared. For this reason it 
is important to characterize BB emission and deter- 
mine what can and what cannot be obtained with it. 
As important as the successes of our model will thus 
be its failures which can guide us toward the correct 
atmosphere, or more generally surface, model (Page 
1995b; Page, Shibanov, & Zavlin 1995). Neverthe- 
less, several of our results are, we hope, sufficiently 
general and robust to be valid despite of the limita- 
tions of BB emission. We restrict ourselves to study 
general properties of the model and refrain from de- 
tailed analyses of the data which we leave for future 
work. Our method to generate surface temperature 
distributions and the observable X-ray fluxes is de- 
scribed in detail in Paper I. 

In § H we present some complements to our sum- 
mary of the ROSAT data given in Paper I. In § || we 
consider the effects of the inclusion of the quadrupolar 
components of the magnetic field and show that when 
they are superposed to a dipolar field they provide a 
sufficiently general configuration for our study. Some 
basic and straightforward results are presented in § ^ 
and ||. We use these configurations to discuss the re- 
liability of our surface temperature model in § ^ and 
also to study the possibility of putting constraints on 
the neutron star (NS) size through gravitational lens- 
ing in § ^| We then discuss the energy dependence of 
the light curves amplitude in § |^, this section will give 
a clear indication of the inadequacy of BB emission 
for understanding detailed features of the observa- 
tions. The next section, ^, adds a second component 
of emission in order to model the higher energy tail of 
Geminga. Section |l^ studies the effect of the surface 
temperature distribution on the boundary condition 
used for modeling the thermal evolution of neutron 
stars. Finally, § [ll] presents our conclusions. 

2. THE ROSAT DATA 

Since the completion of Paper I some information 
has been presented on the second set of ROSAT ob- 
servations of PSR 0656+14 (Ogelman 1995a; Pos- 
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senti, Mereghetti, & Colpi 1996). The complete set of 
data clearly shows the presence of a hard tail which 
is strongly pulsed, as in Geminga and PSR 1055-52. 
Both the energy dependence of the pulsed fraction 
and of the pulse phase show some similarity with PSR 
1055-52 (Ogelman 1995a) and also with Geminga 
(Ogelman 1995b): in all three cases Pf is almost 
constant up to channels ~ 40 (~ 30 for Geminga), 
then decreases and finally increases (very strongly in 
the case of PSR 1055-52). Moreover, at around the 
same energy where Pf increases, the phase of the 
peak changes. The spectral fits with a soft BB for 
the surface thermal emission and a BB or power law 
for the hard tail show that the surface thermal emis- 
sion dominates at channels below channel ~ 50 in the 
case of Geminga (e. g., figure ||) and channel ~ 100 
in the cases of PSR 0656+14 and 1055-52 (see, e. g., 
Ogelman 1995a). An important point is that in the 
case of Geminga the change in the peak phase coin- 
cides with the spectral shift from the surface thermal 
emission to the hard tail (see § ^) and is not surpris- 
ing. However, in the two other cases this phase shift 
apparently, and intriguingly, occurs within the band 
dominated by the surface thermal emission. 

3. QUADRUPOLAR FIELDS 

3.1. The quadrupolar components and grav- 
ity effects 

The first natural step beyond dipolar fields is the 
inclusion of a quadrupolar component which we write 
as 



(i) 



where the five generating fields t>i are listed in ta- 
ble [|. Considering the generating quadrupole compo- 
nents separately at the star's surface, one finds that 
all components reach a maximum strength of 1 and bo 
has a minimum of 1/y/E (=0.45) while the other four 
have a minimum of zero. They have four magnetic 
poles except for bo whose 'south' pole is degenerate 
and becomes a line covering the whole equator. How- 
ever, the general quadrupolar field B Q can have up 
to six poles, south and north poles always being in 
an even number. When the quadrupole is added to a 
dipole we can have again up to six poles, and the num- 
ber of north (south) poles can be odd. Notice finally 
that the scale value of the dipolar component we use 
is the field strength at magnetic pole, i. e., twice the 



value commonly considered, so that all field scales we 
cite always refer to the maximum field strength of the 
cited component at the stellar surface in flat space- 
time. 

We also include the effect of gravity on the mag- 
netic field. It is not as important as red-shift and 
lensing but is nonetheless not completely negligible 
since it can substantially increase the field strength 
and is moreover straightforward to implement. It can 
be written in the form of four multiplicative factors, 
two for the radial component of the magnetic field and 
two for its angular components: if B r , B T (r = 9 or 
4>) denote the radial and angular components of the 
magnetic field in the absence of a gravitational field, 
then: 
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give the dipolar fields in the presence of a gravita- 
tional field, and: 



B% = h 



B Q 



/goo 92 



(3) 



give the quadrupolar fields in the same situation. The 
factors in the previous equalities are given by: 
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where x = Rs/r, Rs = being the star's 

Schwarzschild radius, and <?oo = 1 — a; is the time 
component of the metric (Muslimov & Tsygan 1987). 
These corrections of course tend towards 1 at large 
distances. The values of fi and gi at the stellar sur- 
face, i. e., at x = Rs/R, are plotted in figure [[} the ef- 
fect is stronger on the quadrupole than on the dipole, 
simply due to the stronger r dependence of the for- 
mer. If R — ► Rs then fi and gi tend to oo which 
means that, for a given surface field, the field 'at in- 
finity' vanishes: this is a particular case of the 'no 
hair' theorem for non rotating blackholes (see, e. g., 
Misner, Thorne, & Wheeler 1970). 
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3.2. Statistics of dipole+quadrupole fields 

Since the quadrupolar component introduces five 
more degrees of freedom, a clear general assessment 
of its possible effects is necessary. To investigate this, 
we computed nine sets of 1000 models each, in which 
the quadrupolar components where chosen at random 
by the computer: three orientations of the dipole and 
observer, a = ( = 30°, 60° and 90°, were taken to test 
for geometrical effects (a is the angle between the ro- 
tation axis and the dipole and £ between the rotation 
axis and the observer's direction) and three values for 
the stellar radius, R = 7 km, 10 km and 13 km, were 
used to analyze the effect of gravitational lensing. The 
maximum lensing angles 9 max are respectively, 194°, 
132°, and 117° for radii of 7, 10, and 13 km, respec- 
tively (figure 3 in Paper I). The dipole strength was 
fixed at 10 12 G (value at the magnetic pole in flat 
space-time; GR effects increase it slightly) and the 
1000 quadrupoles were added to the basic dipole. The 
quadrupolc component strengths, Qi, were restricted 
to the range 10 10 G - 5 • 10 12 G. The star's surface 
temperature was calculated using an interior temper- 
ature of T\, = 10 8 K which implies an effective tem- 
perture around 10 6 K (the exact value depending on 
the star size and field configuration) . Notice that GR 
effects depend only on the ratio R/Rs = Rc 2 /2GM 
of the star's radius R to its Schwarzschild radius Rs, 
so that all our results can be extrapolated to other 
masses and radii (see, e. g., figure [l] for easy con- 
version). To save CPU time we only calculated the 
variation of the observable phase dependent effective 
temperature, Tf (Paper I, § 5.2), during the star's ro- 
tation, i. e., the effective temperature of the portion 
of the stellar surface visible to the observer at a given 
phase <fr. Selecting a subset of the 9,000 models so 
generated we performed complete calculations of the 
detectable fluxes to calibrate the observable pulsed 
fraction Pf, in the PSPC's channel range 7 - 50, in 
terms of the pulsation of T*. The results of these 
calculations are shown in figure ^| where each model 
is shown as a black dot and the pure dipole case as 
a black-white ring. The dramatic effect of lensing is 
again clearly seen. For later use, figure ^| plots the 
ratio of the star's effective temperature T e (mag.) for 
the given field configuration to the effective tempera- 
ture it would have without the magnetic field, T e (non 
mag.), at the same interior temperature Tf, = 10 s K. 
A different Tf, would change the values of Pf but not 
change significantly the statistics for the number of 
models giving large Pf. The observable Pf increases 



with decreasing temperature for a given surface mag- 
netic field configuration as explained in § |J. 

We want here to correct our definition of the pulsed 
fraction (Eq. 21 in Paper I) and rather use its stan- 
dard expression as the fraction of counts above mini- 
mum, thus 



(Cts ri 



Cts r 



(8) 



where Cts mean and Cts m i n are, respectively, the 
mean and minimum count rates detected during the 
star's rotation in channel # i. For sinusoidal pulse 
profiles this is equivalent to our former definition and 
does not affect the results of Paper I. 

The first obvious and natural result of figure || 
is that most configurations produce little pulsation: 
most quadrupoles induce several warm regions which 
flatten the pulse profile even below the value of the 
pure dipole case. However, many configurations do 
rise Pf substantially, compared to the purely dipolar 
case, and to a high value in some special cases. We can 
thus obtain observable pulsed fractions higher than 
30% (at energies below 0.5 keV and at T e ~ 10 6 K) 
for a star of radius 13 km. For stars of radii 10 km a 
Pf value comparable to the observed ones can also be 
obtained. The most interesting case is the 7 km radius 
star: despite of the enormous gravitational lensing in 
this star (9 max = 194°), we can still find a configura- 
tion which gives a Pf above 10% at high T e (~ 10 6 K) 
and above 15% at T e ~ 3 • 10 5 K, as we have checked 
explicitly. 

A last point is worth mentioning: the strong in- 
crease in the pulsed fraction possible with the addi- 
tion of a quadrupolar component is not due to an in- 
crease in the area of the cold regions, compared to the 
purely dipolar case, but rather to a displacement of 
the magnetic poles and the surrounding warm regions 
which brings them closer to each other (see, e.g., fig- 
ure ||). The pulse profiles of Geminga and PSR 1055- 
52, which are both considered to be almost orthogo- 
nal rotators, show a single wide peak which indicates 
that the two polar caps and the surrounding warm 
regions are much less than 180° apart. How close the 
magnetic poles can be pushed toward each other is 
what will determine the 'success' of the quadrupolar 
component in increasing the pulsed fraction, whereas 
quadrupoles which moreover introduce several new 
warm regions will fail. Figure |^ illustrates this clearly: 
the ratio T e (mag)/T e (non mag.) is practically the 
same for the strongly pulsed dipole+quadrupole con- 
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figurations as for the purely dipolar case, about 0.9, 
while weakly pulsed configurations which induce sev- 
eral warm regions may have a ratio larger than one. 
These considerations will be of importance for § [l^. 

4. EFFECT OF THE FIELD STRENGTH 

Considering the effect of the field strength, in the 
purely dipolar case we found that the maximum Pf is 
obtained when the overall surface field reaches a few 
times 10 11 G (Paper I): the same result still applies 
to the dipole + quadrupole case. Below 2 — 3 • 10 11 G 
the magnetic effects are getting smaller and above this 
value, even if the magnetic effects are increasing, the 
resulting observable Pf varies very little. The reason 
why the pulsed fraction depends very weakly on the 
field strength for large enough fields is simple: cold 
regions contribute so little to the total flux that it 
does not really matter how cold they are, but only 
how large they are. We are dealing mostly with warm 
plates whose extension is determined more by the field 
orientation than by the field strength as long as the 
field strength does not drop below a few times 10 11 
G. 

5. EFFECT OF THE TEMPERATURE DIS- 
TRIBUTION ON THE SPECTRUM 

The surface temperature distributions induced by 
dipole + quadrupole fields can have quite complicated 
structures. However, as discussed in Paper I for dipo- 
lar fields, the resulting phase integrated spectra are 
close to single temperature spectra at the correspond- 
ing effective temperature. Figure ^ shows examples 
at five temperatures for the dipole+quadrupole field 
from the set generated for figure ^ (see § |3.2j ) which 
gives the highest pulsed fraction. The surface temper- 
ature distribution induced by this field configuration 
is illustrated in figure |. The result is very close to 
the effect of a simple dipole (figure 6 in Paper I) due 
to the fact that the chosen quadrupole induces warm 
regions of approximately the same size as the pure 
dipole but at different locations (see § |3~2| ). When 
taking into account the PSPC's response (figure || B) 
the composite spectrum presents a clear excess com- 
pared to the single temperature spectrum at channels 
above channels 50 - 60. This excess will not affect 
spectral fits in the case of Geminga since the flux at 
channels above 60 is dominated by the hard tail. How- 
ever, in the three other cases where the hard tail only 
appears above 1 keV, even above 1.2 keV for Vela, 



this excess will have some impact on the tempera- 
ture measurement: spectral fits with composite BB 
spectra will give slightly lower r e °°'s than single BB 
spectra and require lower N^'s. (Since the purely 
dipole case gives similar results we must correct the 
corresponding statement made in Paper I about spec- 
tral fits and replace it by the preceding one). We will 
not attempt to quantify this small effect here since 
we restrict ourselves to general characteristics of our 
model. This excess in the Wien tail will also affect 
parameter values for the fits of the hard tail for PSR 
0656+14 and 1055-52 while in the case of Geminga 
our illustrative fit of § ^ can be done with the same 
parameter as used by Halpern & Ruderman (1993) 
for uniform surface temperature. 

6. RELIABILITY OF THE SURFACE TEM- 
PERATURE MODEL 

Our surface temperature model has two ingredi- 
ents: 

1) the value of T s for a given field strength in the two 
cases of parallel and orthogonal transport [T S (Q B — 
0°), and T s (Qb = 90°), respectively, where Q B is the 
angle between the local magnetic field and the radial 
direction] , which are given in figure 1 of Paper I for 
various field strengths B, and 

2) the dependence of T s on Q B for a given field 
strength B: 

T s (e B ) = x (Qb) x T s (e B = o°) (9) 

where 

X (Qb) = (cos 2 e B +xl sin 2 9b) 1 / 4 ; Xo = T °f* = - 

(10) 

As mentioned in Paper I, the dependence given in 
point (2) is a very good approximation, much bet- 
ter than point (1), and its accuracy is certainly bet- 
ter than the statistical uncertainties present in the 
ROSAT data. Most of the uncertainty of the model 
resides in the NS envelope calculations used for point 
(1). The 'T b -T a relationship' T s (T b , B B , B), where T b 
is the interior temperature at the bottom of the enve- 
lope at density p = 10 10 gm cm -3 , is known reason- 
ably well at surface temperatures above 3T0 5 K in the 
case of Q B — 0° and has been calculated by several 
authors who obtained values in reasonable agreement 
with each other; the problem is for the case <d B = 90°, 
where there are large uncertainties due to several fac- 
tors. T S (Q B = 90°) is much lower than T S (S B = 0°), 
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for a given and in this case of orthogonal trans- 
port the envelope structure is often in the "100% 
Murphy regime" where everything that can go wrong 
does go wrong: electron conductivity perpendicular 
to the field is strongly suppressed and photon trans- 
port is the dominant process up to much higher den- 
sities than in the case of parallel transport; at these 
densities and low surface temperatures, the plasma 
temperature is much higher than the local tempera- 
ture but, thanks to the magnetic field, extraordinary 
photons can still exist and transport heat; lattice or 
ion heat transport can also contribute, a phenomenon 
that has been observed in laboratory materials under 
strong magnetic fields (Zyman 1960); the low den- 
sity region (p < 10 4_5 g/cm 3 ) can affect the surface 
temperature and the magnetic effects on the equation 
of state become very large. Finally, meridional flow 
of matter induced by the temperature gradients, ro- 
tation, and magnetic field gradients can also smooth 
the temperature distribution. Unfortunately, none of 
these effects has been taken into account in the en- 
velope calculation that we use for our model (Schaaf 
1990), the only such calculation published to date. 

Despite of these limitations we consider our model 
as quite reliable thanks to the results of § || which 
showed that the observable characteristics of ther- 
mal emission do not depend sensitively on the field 
strength for strengths above a few times 10 11 G 
(when the emission spectra used are from blackbod- 
ies). Maximum magnetic effects are obtained when 
B is superior to about 10 12 G, and saturate at higher 
fields. For a large enough field xo 1 and we can 
thus write 



T S (Q B ) = T s (@b = 0°) x cos 1 / 2 9 2 



(11) 



to a good approximation as long as Os is not too close 
to 90°; when 0# is close to 90° the region is so cold 
that it is not seen anyway. The surface fields of the 
young pulsars we consider can be reasonably assumed 
to be above this threshold value of 2 - 3 • 10 11 G: val- 
ues obtained by the standard magnetic dipole radia- 
tion braking formula, as reported in table 1 of Paper 
I, are between 2.2 • 10 12 G (PSR 1055-52) and 8.8 • 10 12 
G (PSR 0656+14) without GR effects, and are even 
larger when GR effects are included (see above). In 
the case of Geminga the optical data suggest the pres- 
ence of a synchrotron emission line which implies a 
magnetic field of at least 3 • 10 11 G (Bignami et al. 
1996). So, as long as the exact value of T s (Qb = 90°) 
is about 3 times lower than T s (Qb = 0°), as it is at 



3 • 10 11 G in our model, for the field strengths actu- 
ally present, our interpretation of the ROSAT obser- 
vations will remain valid. 

7. GRAVITATIONAL LENSING AND THE 
NS SIZE 

It was shown in Paper I that when only dipolar 
fields are used, the pulsed fraction Pf obtained is 
lower than what has been detected. The addition of a 
quadrupolar component can substantially increase Pf 
to values comparable to the observed ones as shown in 
§ ^2 (figure ||). However, for this to happen we must 
choose carefully the strengths of the various bi com- 
pared to the dipole's strength and orientation: ran- 
dom quadrupolar components generally give tempera- 
ture distributions which are so complicated that they 
flatten the light curves. 

Vela and PSR 1055-52 have pulsed fractions, in the 
channel band where surface thermal emission is de- 
tected, slightly higher than 10% and are considered to 
be almost orthogonal rotators: the models of figure^ 
with a ~ £ = 90° thus apply to them. Taking into 
account that the effective temperature of PSR 1055- 
52 is lower than 10 6 K the observables Pfs would be 
slightly higher than values given in figure f§ We sec 
that for 1.4 M stars with radii of 13 or 10 km many 
models can produce Pfs as high as observed. How- 
ever only very few, special, field configurations can re- 
produce the observed Pf at R = 7 km, less than 0.1% 
of all configurations generated randomly for figure |[ 
The high Pf ~ 20 - 30% of Geminga, also an orthog- 
onal rotator, is however not reproducible with a 7 km 
radius 1.4 M star (i.e., with R = 1.7 Rg): the most 
strongly pulsed dipole+quadrupole model we found 
gives, at T e = 5 • 10 5 K, Pf ~ 15% in channel 
band 7 - 50 as we have checked explicitly, and only 
a few models at 10 km radius reach 20% modulation. 
In the case of PSR 0656+14, which has a BB effective 
temperature of ~ 8 • 10 5 K and a ~ ( ~ 30°, 
we can obtain the observed Pf ~ 14% with a 13 km 
radius 1.4 M star with many models and at 10 km 
radius in only a few cases, but with no model at all for 
a 7 km star: in this latter case the highest Pf we can 
obtain is about 4%. If we take a ~ £ ~ 8° (Lyne & 
Manchester 1988) then Pf is always extremely small: 
even with a radius of 13 km, for a mass of 1.4 Mq , 
the same procedure as used for figure || could not find 
any model which gives a Pf higher that 4%. 
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8. ENERGY AND TEMPERATURE DEPEN- 
DENCE OF THE PULSED FRACTION 

One of the most important characteristics of the 
observed pulsations is the energy dependence of the 
pulsed fraction Pf, We showed in Paper I that BB 
emission with dipolar fields always produces an in- 
crease of Pf with photon energy. We confirm this 
here, with a few exceptions, for more general tem- 
perature distributions and explain the origin of this 
feature. BB emission, with the appropriate surface 
temperature distribution, can explain many of the 
observed characteristics of surface thermal emission 
of the four pulsars we are studying. Since there are 
many theoretical reasons for BB emission to be inade- 
quate it is important to find at least one property not 
reproducible in our simple model in order to guide us 
toward the correct atmosphere model(s). The energy 
dependence of Pf happens to be such a feature: as 
mentioned in § |] in the three cases of PSR 0656+14, 
Geminga, and PSR 1055-52, Pf decreases at channels 
around channel 40 - 50. Another strong constraint on 
atmosphere models is of course that the spectral fit 
must imply a reasonable pulsar distance (cf. discus- 
sion in § p"T| ). 

The reason for the increase of Pf with energy for 
BB emission was elucidated by Page, Shibanov, & 
Zavlin (1995): the hardness of the BB spectrum in- 
creases with energy. The ratio of the BB fluxes emit- 
ted at energy E by two regions of areas A x and A 2 , 
with temperatures T\ and T 2 respectively, is 

F BB {E,T 1 ) = M x exp (E/k B T 2 ) - 1 
F BB (E,T 2 ) A 2 X exp (£/fc s Ti) - 1 [ ' 

which is an increasing function of E if T\ > T 2 . This 
means that if we have a warm spot at temperature 
Ti on a surface with uniform temperature T 2 < T±, 
then Pf naturally increases with energy. A smooth 
temperature distribution with only one peak in tem- 
perature will obviously give the same result. With 
several peaks, by experience we can state that Pf al- 
ways increases (with increasing channel energy) up to 
channels ~ 60 for T e ~ 5 • 10 5 K and up to channel 
~ 100 for T e ~ 10 6 K: we have tried many mod- 
els with dipoles, off-centered dipoles, quadrupoles, 
dipole+quadrupoles, warm plate(s) and found no ex- 
ception. Only in a few special cases, we have 'ob- 
served' a decrease in Pf in the Wien tail which is due 
to the following (figure |J). The surface presents two 
warm regions on opposite sides of the star, a very large 



one and a smaller one where the temperature reaches 
a slightly higher value (figure ||A) , and only one peak 
is visible, i. e., the peak from the small warmer re- 
gion does not appear as a peak but is only filling the 
dip between the successive peaks of the large region 
(figure |^B); by looking at increasing energy the fill- 
ing of the small peak increases (the small region is 
warmer than the large one) and Pf thus decreases 
(figure §C). If the area of the small warmer region 
is increased then it produces a distinctive peak, i. e., 
we obtain a double peak light curve, the amplitude 
of both peaks increases with energy (the amplitude 
of the small region peak will eventually win over the 
large region one) and Pf increases. It seems thus 
that the only mechanism, with BB emission, which 
can produce a decrease in Pf is this filling of the in- 
terpeak by emission from a small warmer region. This 
is in fact the mechanism that Halpern & Ruderman 
(1993) invoked to explain the observed decrease of Pf 
(in channel range 28 - 53 vs. 7 - 28) in Geminga but 
with an essential difference: they proposed that the 
hard tail component, which is about 100° off phase 
with the thermal component, would fill the dip in 
the light curve. However the hard tail is much harder 
than the surface thermal component, i. e., clearly sep- 
arable from it in the spectrum, and cannot actually 
produce a decrease in Pf at energies below 0.5 keV 
as we show in the next section. In the cases where 
we have 'observed' a decrease in Pf, the responsible 
component is only slightly warmer than the rest of the 
star and its emission is practically indistinguishable 
in the spectrum. 

Another feature encountered in Paper I for dipolar 
fields was a steepening of the increase of Pf in the 
channel band 50 - 70 which is due to the response of 
the PSPC. The carbon present in the detector's win- 
dow produces a strong absorption edge at 284 eV: the 
effective area A(E) vanishes just above this energy 
and starts growing significantly only above 400 eV. 
Thus, channels below ~ 40 pick up almost exclusively 
photons with energy below 284 eV while channels 
above ~ 50 detect photons almost exclusively above 
500 eV: since the 'absolute' pulsed fraction, with BB 
emission, is growing with energy there is naturally a 
sharp rise of the observed Pf in the channel band 
50 - 70. This feature must hence be present in any 
model producing a rise of Pf with energy and can be 
seen, for example, in figures |C and @C. 

For a given field configuration, the observable Pf 
at a given energy increases when the overall stellar ef- 
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fective temperature decreases. We stated, wrongly, in 
Paper I that this is due to the increase of the temper- 
ature difference between the warm and cold regions 
when Tb decreases. However, a much stronger temper- 
ature difference increase is induced by the increase of 
the field strength and it has almost no effect on Pf as 
stated in § ||. The actual reason for the increase of Pf, 
at a given energy, with decreasing T e is the decrease 
of Pf with photon energy since lowering T e pushes 
the Rayleigh-Jeans part of the spectrum (which has 
a low Pf) out of the PSPC range. A simple check 
of this statement has been done by considering 'plate 
models' with a uniform surface temperature To on 
which two plates are superposed at temperatures T+ 
and T_, respectively higher and lower than To: keep- 
ing the ratios T + : T : T_ constant but lowering T 
does induce an increase in the observable Pf (at a 
given energy). If we consider the absolute Pf, i.e., 
as would be detected by a detector with absolute en- 
ergy resolution, then the absolute-P/ vs. E curve 
(see for example figure ||C or figure 8 in Paper I) is 
simply shifted leftwards when T is lowered, at con- 
stant T + : T : T_ ratios, and, since Pf increases 
with energy, at a given energy Pf increases when To 
decreases. 

9. SEPARATION OF THE SURFACE THER- 
MAL EMISSION FROM THE HIGH EN- 
ERGY TAIL 

Our study so far, and in Paper I, has concen- 
trated exclusively on the soft X-ray band assuming 
implicitly that this thermal component can be sep- 
arated from the hard tail. We argue here that this 
separation is justified, and show it explicitly in the 
case of Geminga. A simple look at the spectral fits 
(e.g., figure 9A below, or Halpern & Ruderman f993, 
for Geminga; Ogelman f995a for PSRs 0656+14 and 
1055-42) shows that, at energies slightly below the 
crossover of the soft thermal component with the hard 
tail, the contribution of the latter becomes rapidly 
negligible compared to the former. Moreover, the 
spectral fit of the soft component is not changed sig- 
nificantly if the hard tail is fitted by BB emission or 
power law emission (e.g., Halpern & Ruderman 1993). 
One thus does not expect much interference between 
these two components. 

The contribution of the hard tail to the soft band 
is larger if it is modeled as a power law rather than 
a BB. To consider the worst case of interference be- 



tween the two spectral components we thus model 
Gcminga's X-rays with a composite model including 
surface thermal emission and a pulsed power law tail, 
the results being shown in figure ^. The same surface 
temperature model complemented with a hard tail as 
thermal emission from two polar caps has been pre- 
sented in Page & Sarmiento (1996). The spectrum 
(figure ^A) shows that the hard tail emission is about 
30 times weaker than the surface emission at channels 
below 40 but, due to its strong pulsations, its contri- 
bution to the pulse profiles can be about 10% in the 
soft band 7-28 (figure |B1). Nevertheless this is far 
from enough to alter the shape of the light curves in 
both bands 7-28 and 28 - 53 (figure §B1 & B2). 
The important consequence of this is the effect on 
the energy dependence of the pulsed fraction shown 
in figure|]C: the general trend, typical of BB emission, 
consisting of an increase of Pf with energy (see § ||) 
is slightly weakened by the addition of the hard tail 
but not reversed at channels below 100 (figure ||C). 
When modeled as polar cap thermal emission, the 
hard tail contribution to the soft band is almost com- 
pletely negligible (Page & Sarmiento 1996). It is thus 
not possible, within the present framework, to make 
the hard tail responsible for the observed decrease 
of the pulsed fraction with energy below 0.5 keV. In 
conclusion, the observed decrease of Pf with energy, 
at channels below 53, observed in Geminga is most 
certainly an intrinsic property of the surface thermal 
emission which is absolutely irreproducible with BB 
emission. 

The cases of PSR 0656+14 and 1055-52 are more 
delicate but, from a look at the spectral fits there is no 
doubt that the observed decrease of Pf at channels 
below 70 and 50 respectively (Ogelman 1995a) is also 
a property of the surface thermal emission since the 
crossover with the hard tail occurs at channels about 
100 for these two pulsars. However the shift in the 
peak phase (Ogelman 1995a) which occurs within the 
band dominated by the surface thermal emission is 
another feature totally unexplainable with BB emis- 
sion. Thus, both the pulsed fraction variation and the 
peak phase shift must be attributed to the anisotropy 
induced by the magnetic field in the atmosphere, i.e., 
anisotropy of the emitted flux, superposed to a non 
uniform surface temperature distribution. In the case 
of Geminga the shift in the peak phase coincides with 
the shift from the surface thermal emission to the hard 
tail and is not surprising (Ogelman 1995b). 
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10. SURFACE BOUNDARY CONDITION 
FOR NEUTRON STAR COOLING 

A direct application of our surface temperature 
model and of the previous analysis concerns the mod- 
eling of the neutron star thermal evolution. In such 
models, the relationship T b — T s between the tem- 
perature at the bottom of the envelope T b and at the 
stellar surface T s is needed as an outer boundary con- 
dition. Taking into account the surface temperature 
inhomogeneity we must rather speak of a T b — T e rela- 
tionship, where T e is the effective temperature of the 
whole surface, and, given different field structures, 
we can easily generate many such T b — T e relation- 
ships. Cooling models including magnetic effects in 
the envelope have been presented by Van Riper (1991) 
and Haensel & Gnedin (1994) who used a radial field 
T b — T e relationship: this is not a realistic config- 
uration and we will show here that a more careful 
treatment actually leads to conclusions opposite to 
the one reached by these authors. For a radial field, 
T e is increased compared to the non magnetic case, 
at a fixed T b , as shown in figure 1 of Paper I: during 
the photon cooling era, i. e., neutron star age above 
~ 10 5 yrs, this results in an increased surface photon 
emission, for a given internal temperature T b , and ac- 
celerated cooling compared to the non magnetic case. 
For example, with a radial field of strength 10 13 G, 
T e is increased by about 25%, compared to the non 
magnetic case, and thus the photon luminosity is in- 
creased by a factor of 2.5. 

The minimal inclusion of the magnetic field in the 
T b — T e relationship should be done with a dipolar 
geometry. By integrating Eq. |^ for a dipolar field one 
easily obtains 

T e (T b ) = T s (T b , B B = 0) • [1 - 0.47(1 - Xo 4 )] 1/4 (13) 

for the effective temperature T e as a function of the 
maximum surface temperature T s (9s = 0) and xo 
(Eq. |l^): this reduces the magnetic field effects and 
brings the T b — T e relationship close to the non mag- 
netic one as is illustrated in figure fj]. Since we 
have shown that the surface field of the four pul- 
sars which show surface thermal emission (and can 
be compared with cooling models) is not dipolar, we 
must then include the quadrupolar component in the 
Ti, — T e relationship. A general idea of the effect 
can be immediately seen from figure |^ which plots 
T e (mag.)/T e (non mag.). At T b = 10 8 K, i. e., T e ~ 
10 6 K, most configurations also induce a decrease of 



T e compared to the non magnetic case. Moreover, 
T e (mag.)/T e (non mag.) decreases for configurations 
with larger Pf and is systematically smaller than one 
for the cases which can produce the large observed 
pulsed fractions. We also show in figure [?] a second 
T b — T e relationship which corresponds to the most 
strongly pulsed configuration that we have found with 
a dipole+quadrupole field. The interesting result is 
that this relationship does not depend significantly on 
the field configuration: it is very similar to the dipo- 
lar case. We have explicitly verified that other field 
configurations which produce large observable pulsa- 
tions give almost identical results, as can also be seen 
from figure |2j (f or T b = 10 8 K). The discussion at 
the end of § |3.2| anticipated this result and explained 
it beforehand. Configurations of figure || which give 
low P/'s can however increase T e by 10% compared 
to the non magnetic case: the observed strong pul- 
sations, if really due to large surface temperature in- 
homogeneities, totally rule out such field configura- 
tions and thus such boundary conditions. From this 
we conclude that the most likely overall effect of the 
magnetic field, in realistic configurations, is at most 
a slight decrease of the heat flow in the envelope as 
compared to the non magnetic case. Notice that the 
10 13 G case is practically identical to the non mag- 
netic case, a result already obtained by Hernquist 
(1985). Since the strongly pulsed dipole+quadrupole 
field configurations give T b — T e relationships almost 
identical to the purely dipolar case, one could use the 
dipolar relationship for magnetic cooling models. The 
dipolar strengths usually quoted are obtained from 
the standard magnetic dipolar radiation braking for- 
mula and thus only indicative of the actual surface 
field, and even if that formula were exact, GR effects 
can increase the surface field as shown in figure [ij 
We thus propose to take as a boundary condition for 
neutron star cooling models 

T e = a g X s {t T h m , (14) 

where a = 0.85^°;^ and b = 0.55. The value 
of a = 0.85 is an average of the result of Hern- 
quist (1985) who obtained 0.83 and Gundmundsson, 
Pethick, & Epstein (1982) who obtained 0.87 and the 
uncertainties take into account the magnetic field ef- 
fect: +0.05 corresponds to B ~ 10 13 G and -0.20 to 
B ~ 10 11 G. The exponent b = 0.55 is from figure |?] 
and corresponds also to the value of Gundmundsson 
et al. (1982). 

The cooling calculations of Page (1994), restricted 
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to the photon cooling era where the Xj, — T e relation- 
ship is most important, where explicitly performed 
with the non magnetic case in anticipation of the 
present results and his conclusions about the necessity 
of extensive baryon pairing in the core of the Geminga 
neutron star are therefore still valid. For illustration, 
we show in figure || several cooling curves with four 
boundary conditions from figure |^ which show explic- 
itly the smallness of the magnetic field effects. 

When discussing pulse profiles the reliability of the 
envelope models is not very important as discussed in 
§ since most of the effect simply comes from the 
cosr/ 2 &b geometrical factor in equation [H]. When 
considering the T}, — T e relationship the reliability of 
envelope models for Ob = 90° is also not very im- 
portant since the cold regions are so cold that they 
contribute very little to the overall effective tempera- 
ture, [as shown explicitly by the (1 — Xo) dependence 
of T e for dipolar fields in eq. [Uj]. Envelope models 
for Qb = 0° are what really determines the X& — T e 
relationship and more reliable calculations at low tem- 
perature are needed for better interpretation of cool 
neutron star observations which will most certainly be 
detected with the future generations of X-ray satel- 
lites. 

11. DISCUSSION AND CONCLUSIONS 

11.1. Modeling pulsars' surface temperature 
distribution 

We have completed our general study of thermal 
emission from magnetized neutron stars, within the 
framework of blackbody (BB) emission. Our original 
intent was to determine if the inhomogeneous surface 
temperature distribution induced by the anisotropy 
of heat flow in the neutron star envelope is strong 
enough to produce the pulsed fraction observed by 
ROSAT in Vela, Geminga, PSR 0656+14, and PSR 
1055-52. We have shown explicitly in § ^ that most 
of the effect is in the geometrical factor cos 1 / 2 Qb 
(Eq. |ll]), as presented originally by Greenstein & 
Hartke (1983), and fortunately the actual value of xo 
(Eq. |Io|) is of little importance. Moreover, models of 
magnetized neutron star envelopes allow us to relate 
the interior temperature Tb to the maximum surface 
temperature T s (®b = 0°) quite reliably as long as 
the latter is higher than about 3 • 10 5 K. The mini- 
mum surface temperature T s (&b = 90°) is however 
poorly determined but its actual value, equivalent to 
Xo, is not really important. In short, the surface tem- 



perature distribution of a magnetized neutron star can 
be modeled well enough for present practical purposes. 

Since we do not take into account the effect of the 
magnetic field on the emitted spectrum but only on 
the local surface temperature, our results are not sen- 
sitive to the actual strength of the field. As long 
as the magnetic field is stronger than a few times 
10 11 G its effects saturate and we are mostly deal- 
ing with warm plates whose size is determined by 
the field orientation 0#. When atmospheric effects 
are included one can expect a real dependence on 
the field strength due to the presence of absorption 
edge(s) and a field dependent anisotropy of the emis- 
sion. As in the case of dipolar fields (Paper I), the 
composite BB spectrum produced by the temperature 
nonuniformity is close to a BB at the star's effective 
temperature with some excess in the Wien tail (fig- 
ure Q). The quadrupole+dipole configurations which 
induce strong pulsations show an excess similar to the 
pure dipolar case since they have warm and cold re- 
gions of similar areas but many of the configurations 
which produce little pulsations give a composite BB 
spectrum closer to a single temperature BB spectrum 
since they have many warm regions and smaller cold 
regions. 

11.2. The quadrupolar component 

We showed in Paper I that dipolar surface fields 
do not allow to reproduce the amplitude of the ob- 
served pulse profiles; however, we have shown here 
that the inclusion of a quadrupole component is suf- 
ficient to raise the pulsed fraction up to or above the 
observed values. Gravitational lensing is crucial in 
this result, and the effect of the quadrupole, in the 
cases where it does increase the pulsed fraction, is ba- 
sically to push the two magnetic poles, and the warm 
regions surrounding them, closer to each other: the 
areas of the warm and cold regions do not change sub- 
stantially, compared to the purely dipolar case, but 
gravitational lensing is then less effective in keeping 
constantly a warm region in sight. This shift of the 
magnetic poles location is in fact immediately seen 
in the cases of Geminga and PSR 1055-52 which are 
considered to be orthogonal rotators but show only a 
single wide peak in the soft X-ray band while a dipo- 
lar field would produce two peaks (see also discussion 
in Halpern & Ruderman 1993). If we want to stick 
to a multipolar description, strong octupole compo- 
nents would make the matter worse, inducing still 
more warm region and flattening even more the light 
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curves. Thus, the simple observation of significant 
pulsations imply that the surface field is dominated by 
the dipole with a significant quadrupole component, 
but nevertheless weaker than the dipole, and that 
higher order components must be weaker. It has been 
proposed that the presence of the quadrupole is also 
felt in the anomalous dispersion measure of the high 
frequency vs. low frequency radio pulses which has 
been observed in several pulsars (Davies et al. 1984; 
Kuz'min 1992). Similar observations of PSR 0833- 
45, 0656+14, and 1055-52, would allow a direct com- 
parison with the strengths of the quadrupoles needed 
for the interpretation of the radio and X-ray emission 
of these pulsars. No confirmation of the presence of a 
strong quadrupole component could help constraining 
the size of these neutron stars through gravitational 
lensing. Another possibility is that the multipole ex- 
pansion is inappropriate to describe the surface field 
(i. e., would require the inclusion of many strong high 
order components) as would be the case if the field 
has a sunspot like structure: plate tectonics (Ruder- 
man 1991a, b, c) does predict such a structure which 
may be appropriate for cases like Geminga (Page et 
al. 1995). 

Most dipole+quadrupole configurations produce 
little observable pulsations and it is surprising that 
all four observed pulsars have chosen an a priori 
highly improbable configuration which does increase 
Pf. If we visualize the quadrupole as a pair of copla- 
nar dipoles of equal strengths but opposite directions, 
then the dipole+quadrupoles configurations which 
are successful in increasing Pf are almost coplanar: 
there must be a good physical reason for this to hap- 
pen in all four cases. Possible identification of a new 
young neutron star in a supernova remnant close to 
CTB 1 has been recently reported (Hailey & Craig 
1996): this object appears to be similar in age and 
temperature to the Vela pulsar and also shows pul- 
sations at a 10% ± 10% level. If the detected X- 
rays are from the surface thermal emission and Pf 
is above 10% this would raise to five the number of 
neutron stars in which a quadrupole with adequate 
orientation is present. The theory of plate tecton- 
ics (Ruderman 1991a, b, c) predicts that the crust 
of young fast spinning pulsars should break under the 
stress produced by the internal superfluids and super- 
conductors. Plate motion, towards the equator, will 
ensue and the mutual attraction of the two magnetic 
poles may slowly pull them toward each other (as is 
apparently seen in Geminga and PSR 1055-52) induc- 



ing a strong quadrupole component as we need in the 
present study. 

11.3. Constraints on the neutron star size. 

We found that, within our present study, the ob- 
served pulsed fractions of Geminga and PSR 0656+14 
cannot be reproduced if these stars have a radius of 
7 km for a mass of 1.4 Mq , or, more generally, if 
their radius is of the order of 1.7 • Rs, Rs being 
the star's Schwarzschild radius. Moreover, a 10 km 
radius for a 1.4 Mq star, i.e., R — 2.4 Rs, can 
be 'rejected at the 99% confidence level' in the sense 
that less than 10 models out of 1000 generated for 
figure H can reproduce the observed Pf for these 
same two neutron stars. Similarly, for Vela and PSR 
1055-52, R ~ 1.7 Rs can be also be 'rejected at 
the 99% confidence level' for the same reason but 
R = 2.4 Rs is acceptable. Our results thus favor 
rather large radii, above 10 km for a 1.4 Mq mass, 
for all four neutron stars. These conclusions must be 
taken with extreme caution for several reasons: 1) the 
dipole + quadrupole description of the surface field 
may not be adequate, 2) anisotropic emission due to 
the magnetic field may increase significantly the ob- 
servable pulsed fraction, 3) the emitted flux may be 
partially absorbed in the magnetosphere, 4) the sur- 
face temperature may be controlled by factors other 
than heat flow from the hot interior. Each of these 
factors may, however, increase or reduce the observ- 
able pulsed fraction. Thus, we consider our results 
as mostly illustrative but notice that it is encour- 
aging that the observed pulsed fractions fall within 
the range where we could potentially rule out large 
classes of neutron star models through gravitational 
lensing of the surface thermal emission. We finally 
mention that gravitational lensing of the hot polar cap 
emission may also provide complementary estimates 
of the neutron star radius: Yancopoulos, Hamilton, 
k Helfand (1994) find R = (2.4 ± 0.7) R s , or 
10 ± 3 km at 1.4 Mq , for PSR 1929+10 but the 
anisotropy of emission and the lack of information on 
the location of the second polar cap make this esti- 
mate also uncertain. 

11.4. Limitations of blackbody emission 

We confirmed our previous result that with BB 
emission the observable pulsed fraction Pf is un- 
avoidably increasing with photon energy, at low en- 
ergy where the surface thermal emission is detected, 
and explained that this is due to the hardening of 
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the BB spectrum with increasing temperature. The 
three pulsars Geminga, 0656+14, and 1055-52 show a 
decrease of Pf with energy in the soft band, irrepro- 
ducible with BB emission. The model of Page et al. 
(1995) produced such a decrease of Pf by superpos- 
ing two different spectra, a non magnetic hydrogen 
atmosphere at low temperature on most of the stel- 
lar surface and a magnetized hydrogen atmosphere on 
two warm magnetized plates: the magnetic spectrum 
is much softer than the non magnetic one and Pf 
decreases strongly with energy. Extrapolating from 
this result we may speculate that the decrease of Pf 
with energy is due to a softening of the emitted spec- 
trum with increasing temperature. We also showed in 
§ ^| that the separation of the soft thermal component 
and the hard tail is large enough to make the contri- 
bution of the latter in the soft band quite small: the 
observed decrease of Pf is an intrinsic feature of the 
surface thermal emission and not due to an interfer- 
ence of these two components as proposed by Halpern 
& Ruderman (1993) in the case of Geminga. 

A second clear indication for the inadequacy of BB, 
or BB-like, emission is directly found in the spectral 
fits: for a reasonable neutron star size the required 
distance D must agree with other estimates as, e. g., 
the distance obtained from the pulsar's dispersion 
measure or from parallax measurements (Geminga). 
However, the distance D and stellar radius ob- 
tained from spectral fits have to be taken with cau- 
tion due to the uncertainty in the calibration of the 
PSPC in the lower channels band; the fit in this band 
determines the column density Nh which strongly af- 
fects the estimate of the emitted flux and the PSPC 
calibration uncertainty translates into uncertainty on 
D and R°° (Meyer, Pavlov, & Meszaros 1994). In the 
case of Vela, BB emission is disastrous since it implies 
a 3 - 4 km radius neutron star at 500 pc or, alterna- 
tively, a distance of ~ 1500 pc for a standard size neu- 
tron star (Ogelman et al. 1993). The same is appar- 
ently happening in the case of the newly discovered 
neutron star RXJ 0002+6246 (Hailey & Craig 1996) 
for which BB spectral fits also require an effective 
blackbody radius R°° ~ 2.5 — 4.5 km for the esti- 
mated distance of 3 kpc. However, magnetized hydro- 
gen atmosphere models do pass the distance test and 
can fit the Vela spectrum with the standard distance 
of 500 ± 125 pc (Page, Shibanov, & Zavlin 1996) 
and a reasonable hydrogen column density. In the 
case of Geminga, with a parallax measured distance 
of 157±|S pc (Caraveo et al. 1995), available mag- 



netized hydrogen atmosphere models do pretty bad 
since they need a distance of about 20 pc (Meyer et al. 
1994; Page et al. 1995) while BB is acceptable with 
a required distance between 150 and 400 pc (Halpern 
& Ruderman 1993). This low distance needed for 
magnetized hydrogen atmosphere spectral fits may, 
however, be partly due to the uncertainties in the 
PSPC calibration at low energies (Meyer et al. 1994). 
The crucial difference between Vela and Geminga is 
in the temperature: at 10 6 K (Vela) the hydrogen at- 
mosphere is fully ionized and the models are reliable 
while in Geminga, at T ~ 3 — 5 • 10 5 K, the 
atmosphere is only partially ionized and the available 
atmosphere models are not yet very accurate. Atomic 
thermal motion (Pavlov & Meszaros 1993) and its ef- 
fects on the bound-bound (Pavlov & Potekhin 1995) 
and bound-free (Pavlov & Potekhin 1996) opacity will 
affect significantly the absorption edge present at T e 
below 10 6 K, but inclusion of these effects into atmo- 
sphere models has not yet been done and the results 
are hard to predict. The field and temperature depen- 
dence of the absorption edge may naturally induce a 
decrease of Pf with increasing photon energy as ob- 
served (Pavlov 1995). 

These limitations of BB emission may cast doubts 
on the validity of our results. As we stated in the 
introduction we consider that modeling NS thermal 
emission with BB emission is a mandatory first step 
and our study does allow us to delimitate more clearly 
the limitations of this approach. BB is certainly not 
adequate for reliable estimates of pulsars' surface tem- 
perature, as well as size and distance, from spectral 
fits. However, modeling of the pulse profiles is prob- 
ably less dependent on the exact atmosphere charac- 
teristics since the anisotropy of the emission from the 
atmosphere is not very strong at low energy and it is 
unlikely to affect dramatically our conclusions on the 
effect of gravitational lensing and neutron star sizes. 
Our results obviously have to be taken only as a first 
step which hopefully can indicate the direction to be 
taken for future work(s). 

11.5. Neutron star cooling 

We have finally applied our model to assess the ef- 
fect of the crustal magnetic field on the outer bound- 
ary condition used in neutron star cooling models. 
We have shown that the global effect of the magnetic 
field, compared to the non magnetic case, is very small 
and most probably consists in a slight reduction of the 
heat flux in the envelope which results in a lower effec- 
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tive temperature T e for a given inner temperature Tj,. 
This is in sharp distinction to the case where magnetic 
effects are naively applied under the assumption of a 
uniformly magnetized surface, implying an increase of 
T e as compared to the non magnetic case. 

11.6. Warning about the relevance of the 
model 

Finally, we must repeat the warning expressed in 
Paper I about the relevance of our study: strong 
heating of the surface by magnetosphcric hard X- 
rays and/or substantial magnetospheric absorption of 
the emitted flux, as proposed by Halpern & Rud- 
erman (1993), could invalidate seriously our results 
(Page 1995b). The possible detection of pulsed 7- 
rays from PSR 0656+14 (Ramanamurthy 1995) would 
mean that all four pulsars we study are copious 7- 
ray emissors and heating of the neutron star surface 
by such energetic magnetospheres is not an unrea- 
sonable hypothesis. In the original Halpern & Ru- 
dcrman (1993) model for Geminga the flow of elec- 
trons/positrons onto the polar caps produced a lumi- 
nosity L p <~ 2.6 • 10 32 ergs s _1 , radiated away as 
hard X-rays which would then be back scattered onto 
the neutron star surface and reemitted as soft surface 
thermal X-rays. The X-ray luminosity of Geminga is 
about 2 • 10 31 ergs s _1 for the soft thermal com- 
ponent and 8 • 10 29 ergs s _1 for the hard tail (for 
a distance of 160 pc as measured by Caraveo et al. 
1995) so that even a small proportion of the hard po- 
lar X-rays hitting back the surface could explain the 
observed temperature. The problem of this model 
was that the total predicted X-ray luminosity is more 
than one order of magnitude higher than observed. 
The argument has been recently revised and the po- 
lar cap luminosity scaled down by Zhu & Ruderman 
(1996) who now obtain L p ~ 8 • 10 30 ergs s" 1 . 
This is much closer to what is observed and may be 
the factor controlling Geminga's surface temperature. 
However, the soft X-ray flux is about 30 times larger 
than the hard one so that in this model almost all the 
polar cap hard X-rays must be scattered back onto 
the stellar surface. On the other side, neutron star 
cooling models are also perfectly able to explain the 
observed temperature (Page 1994). We can only hope 
that future work(s) will elucidate this dilemma. 
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Table 1 

Spherical components of the five generating quadrupolar fields and location of the magnetic 
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Fig. 1. — Gravitational correction factors for the ra- 
dial (r), fj, and spherical (#,</>), gi, parts of the dipole 
(1 = 1) and quadrupole (1 = 2) components of the 
magnetic field at the stellar surface as a function of 
the ratio of the star's radius R to its Schwarzschild 
radius R$ or of the actual star's radius at masses of 
1.2, 1.4, and 1.6 M . 



Fig. 2. — Statistics of observable pulsations with 
dipolc+quadrupole magnetic fields. For each stellar 
radius and dipole-observer orientation (a-£) with re- 
spect to the rotational axis 1000 surface magnetic field 
and the induced temperature distribution were gen- 
erated The x axis plots the observable pulsed fraction 
in the PSPC channel band 7-50 for the given sur- 
face temperature distribution and the y axis the ratio 
of the star's effective temperature for the given mag- 
netic field configuration to the temperature it would 
have with zero field. The surface temperature are cal- 
culated from an interior temperature Tf, = 10 8 K: all 
models have an effective temperature around 10 6 K. 
Lowering T e does increase Pf: by about 15% at 
T e ~ 3 • 10 5 K and 30% at 10 5 K. See text, § |J, 
for more detail. 



Fig. 3. — Surface temperature distribution for the 
dipole + quadrupole magnetic field configuration of 
figure H which gives the highest observable pulsed 
fraction. The whole star's surface is shown in an area 
preserving mapping. With an internal temperature 
of 10 8 K, the effective temperature is 9.33 • 10 5 K 
while the maximum and minimum temperatures are, 
respectively, 1.2 • 10 6 K and 3.1 • 10 5 K (1.4 M star 
with a 10 km radius). Four sets of isotherms are plot- 
ted as indicated. Field components: dipolar compo- 
nent of strength 10 12 G and orientation 9 = 90° and 
= 0°, and quadrupolar components of strengths 
Qo = 2.48 • 10 10 G, Qi = 6.84 • 10 10 G, 
Q 2 = 5.01 • 10 10 G, Q 3 = 4.80 ■ 10 11 G, 
and Q4 = — 6.24 ■ 10 10 G: the maximum field 
strength at the surface is 2.18 • 10 12 G and the min- 
imum 1.72 • 10 11 G. The two black dots show the 
approximate locations of the polar caps, obtained by 
integrating the last open magnetic field lines from the 
light cylinder down to the stellar surface while the two 
small circles at 9 = 90° and = 180° and 360° show 
their positions without the quadrupolar component. 



Fig. 4. — Effect of the inhomogeneous surface tem- 
perature on the spectrum. A: incident flux (i.e. with 
redshift and lensing) at five temperatures, with zero 
interstellar absorption. B: the same fluxes with in- 
terstellar absorption and passed through the PSPC 
response matrix. The temperatures of 5, 7, 9 and 15 
•10 5 K correspond to the estimated blackbody tem- 
peratures of Geminga, PSR 1055-52, PSR 0656+14 
and Vela respectively. We used a 1.4 M Q , 10 km ra- 
dius (R°° = 13.06 km) star, at a distance of 500 pc, 
with the same field configuration as in figure |^. Con- 
tinuous curves: composite BB spectra, dashed curves: 
single temperature BB spectra at T = T e . 



Fig. 5. — A model producing a decreasing pulsed frac- 
tion with blackbody emission. (1.4 M Q star with a 12 
km radius.) 

A. Surface temperature distribution. Magnetic field: 
dipole orthogonal to the rotation axis (in direction 
9 = 90° and = 90°) with its center shifted by a half 
radius in direction of the north magnetic pole. The 
warm region around the north pole reaches a temper- 
ature of 8.15 • 10 5 K but is much smaller than the 
warm region around the south pole whose maximum 
temperature is 7.25 • 10 5 K. 

B. Light curves in three channel bands produced by 
the temperature distribution of A, as seen through 
the ROSATs PSPC by an observer at 150 pc and 
Nh — 10 20 cm~ 2 . Star's rotation axis along 9 = 0° 
and observer at 9 = 90°. The dotted lines show light 
curves where the temperature of the smaller north 
magnetic pole region has been blocked at 7.25 • 10 5 
K, i. e., the same maximum temperature as the large 
region around the south magnetic pole: this clearly 
shows that this warmest region is responsible for the 
decrease of the Pf. 

C. Observable (with the PSPC) and absolute (i. e., as 
would be seen with a perfect energy resolution detec- 
tor) Pf. The dotted curves show the corresponding 
Pf when the temperature of the small warm region 
has been blocked at 7.25 • 10 5 K (as in B). 
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Fig. 6. — Composite model of Geminga's soft X- 
ray emission: surface thermal emission (dashed lines) 
plus power law emission (dotted lines) and total emis- 
sion (continuous lines). Surface temperature dis- 
tribution induced by a dipole+quadrupole magnetic 
field: dipole of strength 10 12 G perpendicular to 
the rotational axis and quadrupole with components 
Qo = -2.28 • 10 11 G, Qi = +4.95 • 10 10 G, 
Q 2 = +3.36 ■ 10 10 G, Q 3 = +7.08 ■ 10 11 G, 
Q4 — — 4.53 ■ 10 11 G. The interior temperature 
is 4.04 • 10 7 K, giving an effective temperature at 
infinity T e °° = 4.30 • 10 5 K for a 1.4 M star with a 
radius of 12 km. The hard tail is a sum of two compo- 
nents, both power laws with indices 1.47, modulated 
with the rotational phase $ by a cos 4> factor. The ob- 
server is located within the rotational equatorial plane 
at a distance of 185 pc and N H = 2.1 • 10 20 cm~ 2 . 
A. Observable count rates from the surface thermal 
emission, the power law emission, and the total emis- 
sion compared with the ROSAT data (from Halpern 
& Ruderman 1993). B. Pulse profiles in three chan- 
nel bands, compared with the ROSAT data (from 
Halpern & Ruderman 1993). C. Pulsed fraction for 
the two components and the total emission. 



Fig. 7. — Tf, — T e relationship: thick lines correspond 
to purely dipolar fields while the thin lines correspond 
to one of the maximally pulsed dipole+quadrupole 
configuration of figure ||. The indicated field strengths 
are the dipole component strength at the magnetic 
pole (without GR correction: so the actual field value 
is about 50 - 100% higher, see figure pi). 



Fig. 8. — Effect of the crustal magnetic field on the 
cooling of neutron stars. We show the 'standard' 
cooling of a 1.4 M Q neutron star built with the FP 
equation of state (Friedman & Pandharipande 1981) 
with and without core neutron 3 P2 pairing. The neu- 
tron core pairing critical temperature T c is taken from 
Hoffberg et al. (1970): the high T c of this calculation 
shows the maximum effect of pairing. Crust neu- 
trons are always assumed to be paired in 1 Sq state 
with T c taken from Ainsworth, Wambach & Pines 
(1989). Core neutron pairing has no effect at the 
early stage (age below ~ 20 yrs) but (i) slows down 
the cooling during the neutrino-cooling era (age up to 
~ 10 5 — 10 6 yrs) by strongly suppressing the neutrino 
emission, and (ii) fastens the cooling later during the 
photon-cooling era by suppressing the star's specific 
heat by about 75% (see, e. g., Page 1994). The two 
sets of four curves show the effect of the crustal mag- 
netic field on the cooling by its influence on the outer 
boundary condition. The four T5 — T e relationships 
are taken from figure ^. During the neutrino-cooling 
era the surface temperature follows the evolution of 
the core temperature: a Tf, — T e relationship which 
gives a higher T e , for a given Tf,, results naturally in 
a higher stellar effective temperature. During the fol- 
lowing photon-cooling era the effect is inverted since 
the dominating energy loss is from surface thermal 
emission with a luminosity L 1 = ^nR 2 aT^ cx T 22 
(see equation [Ii]) and a higher T e for a given XJ, re- 
sults in a larger L™. Calculations were performed with 
the cooling code used by Page (1994) 
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